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Abstract 

We search for the doubly charmed baryonic decay B° — > A+A~, in a data sample of 520 x 10 6 
BB events accumulated at the T(45) resonance with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider. We find no significant signal and set an upper limit of B(B° — > A+A~) < 
6.2 x 10~ 5 at 90% confidence level. The result is significantly below a naive extrapolation from 
B(B~ — ► H^A^) assuming a simple Cabibbo-suppression factor of \V c d/V cs \ 2 . The small branching 
fraction could be attributed to a suppression due to the large momentum of the baryonic decay 
products, which has been observed in other charmed baryonic two-body B decays. 

PACS numbers: 13.20.He 
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The large mass of the b quark and the large quark mixing matrix element V c b [H, 0] for 
the b — > c transition lead to a large branching fraction (~ 10%) 0] for charmed baryonic 
decays of the B meson. Charmed baryonic decays into four-, three- and two-body final 
states have already been observed. The measured branching fractions; B(B° — > A + p7i + 7i~) = 
(1.12±0.05±0.14±0.29) x 10" 3 fl, B{B~ -> A + pn~) = (2.01 ±0.15±0.20±0.52)_x 10" 4 Q 
and B(B° -> A + p) = (2.19±g;^ ± 0.32 ± 0.57) x 10~ 5 [6] (also see Ref. @, B I, Ej), point 
to a hierarchy of branching fractions depending on the multiplicity in the final state [HI ]. 
The measurements provide stringent constraints on theoretical models for charmed baryonic 



decays of the B meson |12|, [13j, 114 



The hierarchy can be understood by large contributions of various intermediate states 
known in the decays j3, 0, @, 0, 10]. The key is to understand quantitatively the decay 
mechanism of the two-body decays. For example, B(B~ — > S c (2455)°p) = (3.7 ± 0.7 ± 0.4 ± 
1.0) x 10~ 5 [5j is observed in the three-body decay B~ — > A + pir~, which is comparable to 
B(B° — > A + p ). There is also an interesting measurement of B(B~ — > 5°A~) = (5.8 ± 2.3) x 
10~ 3 [l5|, 0, [Dl H, 3], which is quite large in comparison with B(B° — ► Aj 



\p), and does 

not follow the hierarchy. Figures QJa) and (b) show quark diagrams relevant for these decays 
through Cabibbo-favored b — > cW~ transitions with W~ — > ud and W~ — > cs, respectively. 
Since we naively expect similar branching fractions as IV^V^I 2 ~ lKt^cs| 2 , the two-order of 
magnitude difference between B(B~ — > S°A~) and B(B° — > A+p) is a puzzle. It indicates 
that there is some mechanism to enhance or suppress specific two-body decays. A discussion 
of a dynamical suppression mechanism, based on the large Q-value in B° — > A+p compared 
to B~ — > S^A^, is given in Ref. [2(|. It is important to study various two-body decays to 
understand charmed baryonic B decays. 




a ) A oc V* b V ud 



b) a oc vsy cs 



c) A oc VIVcd 



FIG. 1: Quark diagrams for (a) B° -» A+p, (b) B~ -* E° C A~ and (c) B° -> A+A~. The first two 
decays are Cabibbo-favored with CKM couplings V* b V u d and V^V CS , respectively, while the third 
one is Cabibbo-suppressed with coupling V* b V C( i- 



In this report, we study the doubly charmed baryonic decay B° — > A+A~ as shown 
in Fig. die). This mode is naively expected to have a branching fraction suppressed by a 
Cabibbo factor of 5.4% [3| relative to the Cabibbo-favored decays. Given the large branching 
fraction B(B~ — > H^A^) relative to B(B° — > A+p), we search for the decay B° — > A+A~ 
and compare the observed branching fraction with simple estimates. We expect B(B° — > 
A+A ( 7) = (7.7±3.0) x 10 -7 from B(B° — > A+p), taking into account the Cabibbo-suppression 
factor and the phase space factors in two-body decays proportional to the decay momentum 
in the B rest frame (assuming a relative S-wave, L=0). Alternatively, we expect B(B° — > 
A+A~) = (3.6± 1.1) x 10~ 4 from B{B~ — > 2°A^r). We expect 0.1 and 45 events, respectively, 
from these two estimates scaled to our data sample 
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This analysis is based on a data sample of 479 fb _1 , corresponding to 520 x 10 6 BB 
events, which were recorded at the T(4S) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider [2l| . 

The Belle detector is a large-solid- angle spectrometer based on a 1.5 Tesla superconduct- 
ing solenoid magnet. It consists of a three layer silicon vertex detector for the first sample of 
152 x 10 6 BB pairs, a four layer silicon vertex detector for the later 368 x 10 6 BB pairs, a 50 
layer central drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), 
a barrel-like arrangement of time of flight scintillation counters (TOF), and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals located inside the superconducting solenoid 
coil. An iron flux return located outside the coil is instrumented to detect K\ mesons and 
to identify muons. The detector is described in detail elsewhere (22|. To simulate detector 
response and to estimate efficiency for signal measurement, we use Monte Carlo (MC) event 



generation program EvtGen [23j and a GEANT [24j based detector simulation code. 

To search for B° — > A+A~ we reconstruct a pair of A+'s decaying into pK~n + . Charge- 
conjugate modes are implicitly included throughout this paper unless noted otherwise. We 
require tracks to have a distance of closest approach to the interaction point less than 5 cm 
along the z-axis (opposite to the e + beam direction) and 1 cm in a plane perpendicular to the 
z-axis. Hadrons (protons, kaons and pions) are identified by using likelihood ratios based 
on CDC dE/dx, TOF and ACC information. We use likelihood ratios L S /(L S + L^), where 
s and b stand for the hadron species to be identified and for the others, respectively. We 
require the ratios to be greater than 0.6, 0.6 and 0.4 for proton, kaon and pion selection, 
respectively. The efficiency for proton identification is 95% with a kaon fake rate of 1.0% 
due to the small proton momentum (~ 1 GeV/c) in these baryonic decays. The efficiencies 
for kaons and pions are about 90%, while the corresponding pion and kaon misidentification 
rates are approximately 10% |25J. Tracks that are positively identified as electrons or muons 
are rejected. We impose loose requirements on the vertex fit x 2 's for A+ — > pK~ir + (x^+) 

and B° — > A+A~ (xb) to reject background from the decay products of Kg and A particles. 
When there are multiple B candidates (3%) in an event, we choose the candidate with the 
smallest \\. 




FIG. 2: A+(pK~-K + ) mass distribution for B° — >■ A+A~ candidates in \AE\ < 0.2 GeV and 
5.2GeV/c 2 < M bc < 5.3GeV/c 2 . (a) Data and (b) MC signal. The curves show the fits with a 
double Gaussian for the signal and a linear function for the background. 
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We search for the B signal in the two dimensional plane of AE and M bc . The variable 
AE = Eb — -Ebcam is the difference between the reconstructed B meson energy (Eb) and 



the beam energy (-E bea m)- ilf bc — \J E^ e&m — P B is the beam energy constrained B meson 
mass with the momentum vector of the B meson (Pb)- Here -E bea m, Eb and Pb are defined 
in the center-of-mass system (CMS). We use the A+ mass [s| and the measured momentum 
of the A+ system to calculate Eb, as it gives a better AE resolution, 4.3MeV/c 2 , than 
that calculated with the A+ energies reconstructed from the decay products, 6.6MeV/c 2 . 
To optimize the selection parameters for the signal search, we define a B signal region of 
\AE\ < 0.02 GeV(~ 4a) and 5.27 GeV/c 2 < M bc < 5.3 GeV/c 2 . 

Figure [2] shows the A+ mass distribution for (a) data and (b) the MC signal for B signal 
candidates with \AE\ < 0.2 GeV and 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 . We find a significant 
A+ mass peak in the data due to the large inclusive branching fraction for B meson decays 
with a A+ baryon in the final state. The curves show fits using a double Gaussian for the 
signal and a linear function for the background. We obtain a A+ yield of 1281 ± 69 events 
with a x 2 /ndf = 59.4/65 (67.4%). In the fit to the data, we fix the ratio of cxtaii/ccorc to 2.29 
and the tail fraction (to the total area) to 0.284; these values are obtained from a fit to the 
MC signal. The parameters a ta u and cr cor e are the widths for the core and tail Gaussians, 
respectively. The fitted masses and a core are (2285.3 ± 0.2) MeV/c 2 and (3.3 ± 0.2) MeV/c 2 
for the data, and (2285.9 ± 0.1)MeV/c 2 and (3.2 ± 0.1)MeV/c 2 for the MC signal. We 
require that the A+ masses lie in the range 2.275 GeV/c 2 to 2.295 GeV/c 2 (±3a core ). The 
small differences between the data and MC signal are taken into account in the systematic 
error as discussed later. 

In this analysis, the A+ mass requirements are very effective in suppressing the continuum 
background (e + e~ — > qq,q = u, d, s, c). The dominant background is from generic B events. 
To suppress the background further, we use the variable cos#b, which is the cosine of the 
angle between the reconstructed B direction and the e~ beam direction in the CMS. The B 
signal has a (1 — cos 2 #b) distribution while the generic B background and the continuum 
background have a nearly flat distribution. Using MC simulation, we examine the figure 
of merit S/VS + N function of cos#b- Here, S and N are the signal and background 
yields in the B signal region, respectively. We assume a branching fraction B(B° — > A+A ( T) = 
5 x 10~ 5 and a sample of 6 x 10 8 BB events, and optimize the figure of merit with the 
requirement | cos^b] < 0.8. 

To obtain the signal yield, we perform an unbinned maximum likelihood fit to the 
B° — > A+A~ candidates in a two-dimensional (2D) region —0.15 GeV < AE < 0.2 GeV 
and 5.2GeV/c 2 < M bc < 5.3 GeV/c 2 . We exclude the region AE < -0.15 GeV, as we 
find from MC simulation that a background from B~/° — > A+A~7r~/° populates the region 
AE ~ —0.2 GeV. Thus, the effect of the background is negligibly small (<0.05 events) in 
the fit region, even if we assume large values of B(B~/° — > A+A~ K~/°) 

We use a likelihood defined by 



e -(n s +n b ) _n 



n - i=l 



Y[[n s F s (AE h M hci ) + n h F h (AE h M hci )} (1) 



with the signal yield n s and the background yield n b . The parameter n is the observed 
number of events. The probability density function (PDF) for the signal F s (AE,Mb c ) is 
expressed as a product of a double Gaussian in AE and a single Gaussian in M bc , while the 
PDF for the background F^(AE, M bc ) is expressed as a product of a linear function in AE 



and an ARGUS function 27f in M bc 



6 



In the fit, the AE and M bc signal shape parameters are fixed to those obtained from 
one-dimensional fits to the individual simulated distributions for AE with 5.27GeV/c 2 < 
M bc < 5.30GeV/c 2 , and M bc with \AE\ < 0.02 GeV. The yields n s and n h , the AE linear 
slope parameter and the ARGUS shape parameter are floated. We obtain a signal efficiency 
of 0.106 ± 0.001 from a 2D fit to the MC signal. For the fit to the data, we fix the signal 
parameters to those calibrated for the MC/data systematic difference by using a control 
sample of B° — > A+]97r + 7r~ decays. 




FIG. 3: Two-dimensional unbinned likelihood fit to the data in -0.15 GeV < AE < 0.20 GeV 
and 5.20GeV/c 2 < M bc < 5.30GeV/c 2 . (a) AE and (b) M bc distributions for all events, (c) AE 
distribution for 5.27GeV/c 2 < M bc < 5.30GeV/c 2 , and (d) M bc distribution for \ AE\ < 0.02 GeV. 
(e) AE distribution for M bc < 5.27GeV/c 2 and (f) M bc distribution for \AE\ > 0.02 GeV. The 
curves represent the fitted signal (dotted lines) and the total (solid lines) yield. 

Figure [3] shows the fit to the data. We obtain a signal of 2.7i 2 ; events with a statistical 
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significance of 1.6a. The significance is calculated as y — 21n(L /L max ), where L max and L 
are the likelihood values at the fitted signal yield and the signal fixed to zero. 

We investigate a possible peaking background in the sideband data, which includes a 
background from B° — > A~^pir + 7C~ [4], when a tx + is misidentified as a K + . We define the 
sideband by requiring that one of the A+ candidate masses lies in the range 2.245 GeV/c 2 — 
2.325 GeV/c 2 while excluding masses in the range 2.275 GeV/c 2 — 2.295 GeV/c 2 . From the 
2D fit to the sideband, we estimate a peaking background of —0.1 ± 0.5 events, which is 
consistent with zero. 

We estimate a systematic error of 14.5% in event reconstruction and selection; a 12.6% 
uncertainty in the efficiency (arising from possible differences between the data and MC 
simulation in the reconstructed A+ mass, particle identification and tracking), a 7.1% un- 
certainty due to the uncertainty of the signal parameterization used in the 2D fit (obtained 
by varying the parameters by one standard deviation), and a 1.3% uncertainty in the total 
number of BB events. We obtain a total systematic error of 62% in the measured branching 
fraction, including a 58% uncertainty due to an error in £>(A+ — * pK~ix + ) = (5.0 ± 1.3)% |3j 
and a 18% error for the peaking background. We correct the signal efficiency by a factor 
of 0.90 due to a systematic difference in particle identification between MC and data. We 
assume the same numbers of neutral and charged BB pairs, and obtain a branching fraction 
of (2.2j£g(stat) ± 1.3(syst)) x 10" 5 . 

We calculate 7.7 events for the upper limit yield at 90% confidence level (CL) by 
integration of the likelihood function obtained from the 2D fit. We use the formula 
of 90% = Jq ul L(n|s)cfe// °° L(n\s)ds with n = 2.7, where the likelihood L(n\s) = 
ff^ L fit (n|s*) • G(s — s*)ds* is convolved with the Gaussian G(s — s*) to take into account 
the total error, which is composed of errors in the fitted yield (the signal and the peaking 
background), and the systematic error discussed above. The corresponding upper limit is 
found to be B(B° -* A+A") < 6.2 x 10~ 5 at 90% CL. 

The present result is much smaller than the naive estimate of (3.6 ± 1.1) x 10~ 4 from 
B(B~ — > S°A~) [l5[ with a significance of approximately 3a, where the main uncertainty 
comes from the experimental error in B(B~ — > S°A C 7)-S(H° — > H~7r + ). On the other hand, 
no significant difference is observed for the naive estimate of (7.7±3.0) x 10 -7 from Br(B° — > 
A^p) |6j due to the limited statistics. Figure H] compares the result with the data for other 
charmed baryonic two-body B decays; B~ — > 2°A~, B~ — > S c (2455)°p [5] and B° — > A+p. 
We define a rescaled branching fraction T = B/ (p-CSF) [23]. Here p is the decay momentum 
in the B rest frame, which represents a phase space factor (assuming L=0), and CSF is 
a Cabibbo suppression factor [3]: 1.0 for B~ — > H°A~ and B° — > A^p, and 0.054 for 
B° -> A+A~. We also plot T{pp) UL for the 90% CL upper limit on B(B° -> pp) Q with 
CSF=|K 1 fc/V^6| 2 = 0.011 [3] assuming a b — ► u{du) tree transition. The open and solid points 
with error bars show the data for B~ and B° decays, respectively. The dashed line shows 
the function ln(jF(p)) = c + s x p with s = —6.9 ± 0.8 (GeV/c)" 1 to guide the eye, which is 
obtained by a fit to the three data points. The 90% CL upper limit ^(A+A")^ is close to 
the line. 

In summary, we search for the doubly charmed baryonic decay B° — > A+A~ in a data 
sample of 520 x 10 6 BB events. We obtain B(B° -> A+A") = (2.2l 2 J(stat) ± 1.3(syst)) x 10~ 5 
with an upper limit of B(B° — > A+A^) < 6.2 x 10~ 5 at 90% confidence level. The result 
is significantly smaller than a naive extrapolation from B(B~ — > H^A^), assuming a simple 
Cabibbo suppression factor. The suppression of B° — > A+A~ could be attributed to the 
strong momentum dependence of the decay amplitude that has been observed in other 
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FIG. 4: The rescaled branching fraction T = B/(p • CSF) for B~ -» H°A~, B° -> A+A~, 
B~ ^> E^p, 5° — > A+p and S° — ► decays. The dashed line shows a fit to ln(.F(p)) = c|sxp 
with s = —6.9 ± 0.8 (GeV/c) _1 to guide the eye. 



charmed baryonic two-body B decays. 
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